Effect of magnetic frustration on single-hole spectral function in the t-t'-t"-J model 
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We examine the effect of the magnetic frustration J' on the single-hole spectral function in the 
t-t-t -J model. At zero temperature, the exact diagonalization (ED) and the self-consistent Born 
approximation (SCBA) methods are used. We find that the frustration suppresses the quasiparticle 
(QP) weight at small momentum k, whereas the QP peak at k=(7r/2,7r/2) remains sharp. We also 
show the temperature dependence of the single-hole spectral function by using the ED method. It is 
found that the lineshapes at (n/2,0) and (7r/2,7r/2) show different temperature dependence. These 
findings are consistent with the angle-resolved photoemission data on Sr2Cu02Cl2, and indicate the 
importance of the magnetic frustration on the electronic states of the insulating cuprates. 

PACS numbers: 75.50.Ee, 78.20.Bh, 78.66.Nk, 79.60.-i 



I. INTRODUCTION 

Since the discovery of high temperature superconduc- 
tivity, the t-J model has been extensively studied as 
a theoretical model for the superconductors. However, 
angle-resolved photoemission spectroscopy (ARPES) 
measurements on Sr2Cu02Cl2 by Wells et alu have 
clearly shown that the t-J model does not explain the 
experimental data near k=(-7r,0): the t-J model predicts 
that the energies of the quasiparticle (QP) at (ir,0) and 
(7r/2,7r/2) are very similar, while the experimental energy 
at (7r,0) isJowpr than that at (ir/2,ir/2). Recent theoret- 
ical workstatLil have shown that this discrepancy can be 
resolved by introducing hopping matrix elements to sec- 
ond and third nearest neighbors (t 1 and t") into the t-J 
model. The Hamiltonian, termed the t-t'-t"-J model, is 
given by 

Htt't" j = Htt't" + Hj , (1) 
Htt't" =-* J] c\ a cja -t' S L S i° 

-*" E 4^+H.c. , (2) 
Hj = J J2 S i- S J . ( 3 ) 

(id) 1st 

where ci a = Ci CT (l — ni_ CT ) is the annihilation operator 
of an electron with spin a at site i with the constraint 
of no double occupancy, Si is the spin operator and the 
summations (i,j) lst , ano - (*>i)3rd run over nrs t, 

second and third nearest-neighborj-pairs, respectively. 

Very recently, it has been argueda that not only the QP 
dispersion but also the lineshape of the ARPES data con- 
tains valuable information about the electronic states. It 
was shown that the t-t'-t"-J model describes the observed 
broad peak at (ir,0). A weakening of the antiferromag- 
netic (AF) spin correlation induced by the charge-carrier 
motion due to the longe-range hoppings is found to be 
responsible for the broad lineshape. 



The fact that the lineshape of ARPES spectra is sen- 
sitive to the AF correlation in the spin background sug- 
gests that, if magnetic frustration in the spin system is 
strong, the lineshape of the spectra should be also af- 
fected significantly. A dominant interaction which causes 
the magnetic frustration is considered to be the next- 
nearest neighbor exchange interaction J' . The Hamilto- 
nian is given by 

Hj, = J' ]T 3i-3 s . (4) 

(id) and 

The two-dimensional (2D) spin-| Heisenberg Hamilto- 
nian including both J and J', i.e. Hj + Hj>, has, been 
extensively studied, motivated by the suggestionE£l that 
the physics of the high-T c phenomena is closely related to 
the presence of non-Neel states in insulating 2D cuprates. 
The phase diagram of the-J^J' model obtained as a func- 
tion of J'/ J is as follows £3'Ej the Neel state is the ground 
state for 0< J'/ J<0.6, while a collinear state becomes 
stable in the ground state for 0.6<J'/J<1. Since the 
realistic value of J' j J in ihe 2D insulating cuprates is 
expected to be 0.1~0.2,ll2~E3 the ground state remains to 
be the Neel state. Although the Neel state is stable, the 
excited states obtained after kicking photoelectron out of 
the system might be sensitive to the magnetic frustration. 
As far as we know there has been no study to examine 
the effect of J' on the single-hole spectral function. 

In this paper, we examine the effect of the magnetic 
frustration J' on the single-hole spectral function at 
zero and finite temperatures in the t-t'-t"-J model. At 
zero temperature, the exact diagonalization (ED) and 
self-consistent Born approximation (SCBA) methods are 
used. The frustration is found to suppress the QP weight 
at small momentum. The suppression of the QP peak 
gives rise to a broad spectrum at the momentum, whereas 
the QP peak at (ir/2,ir/2) remains sharp. This is con- 
sistent with the ARPES data on S^CuC^C^ that the 
(tt/2,0) spectrum is broader than the (tt/2,tt/2) one. We 
also calculate the single-hole spectral function at finite 
temperatures in the t-t'-tf'-J model with J' by using the 
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ED method. It is shown that the lineshapes between the 
(n/2,0) and (w/2,ir/2) spectra behave differently with in- 
creasing temperature. This difference is also consistent 
with experimental data by Kim and ShenEj These find- 
ings indicate the importance of the magnetic frustration 
on the electronic states of the 2D insulating cuprates. 

II. METHOD OF CALCULATIONS 



and l^>r]>0. The AF long-range order is assumed in the 
derivation of these equations. In the next section, wc 
will compare the SCBA results with ED ones and see 
the validity of the assumption. A (k, to) is assumed to be 
equal to the imaginary part of the hole Green's function 
Eq. (7): 

A(k,uj) = —TmG(k,u) . (15) 



The single-hole spectral function at temperature T is 
given by 

4(k, W ) = I^|(^- 1 |c kCT K>| 2 



x5(u>-E%+EZ-*)e 



-PK 



(5) 



where 4^ is the wave function of the m-th eigen- 
states with energy E% in the iV-electron system. Z = 

J2 m e ~ f3E ™ 1S tne partition function with f3 — l/k^T. 
The Boltzmann factor fcs is taken to be 1, hereafter. The 
initial states are the eigenstates of the J- J' Heisen- 
berg model. At T—0, the initial state is restricted to 
the ground state of the J- J' model. In the evaluation of 
Eq. (g) at finite temperatures, we use a 4x4 lattice with 
periodic boundary condition. After diagonalizing the J- 
J' Hamiltonian by the Householder method, we make use 
of the Lanczos procedure for each initial eigenstate in the 
Hamiltonian to obtain A(k,u>). 

At T=0, we also use the self-consistent Boriii-approx- 
imation (SCBA) method to calculate ^4(k, w).ETlj In 
the SCBA method, the hole Green's function G (k, ui) 
and the self-energy S (k, iS) satisfy a set of self-consistent 
equation for a 2D square lattice with N sites: 
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G(k, W ) = 

LU — L (k, LO) — £k + tT) 

E(k,w) =^F(k,q)G(k-q,c^-^ q ) , 
q 

32t 2 

-^(k,q) = |7k-qUq + 7kW q | 2 , 

where 

ek = 4t' cos k x cos k y + 2t" (cos 2k x + cos 2k y ) 



Wq = 2Jyjq 2 - 7q 2 , 
J' , 

j q = 1 - — (1 - cosq x cosq v ) , 
7k = 2 {cosk x + cosk y ) , 



V2 \ V M 2 - 7q 2 



-sign (7 q ) -= 



V2 \ V V - 7q 2 



(6) 
(7) 

(8) 

(9) 
(10) 

(11) 
(12) 

(13) 
(14) 
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FIG. 1. Single-hole spectral function A(k, cj) of the 
t-t'-t"-J model with various values of J' at T=0. t'=-0Mt, 
i"=0.23t and J=0At. J' /J is (a), 0.1 (b), 0.3 (c), and 0.5 
(d). Calculations are performed on the 4x4 cluster by the ED 
method (solid line) and the SCBA method (dashed line). The 
SCBA data at (n,n) and (n,n/2) are excluded from the pan- 
els, since they are equivalent to those at (0,0) and (n/2,0), 
respectively. The solid lines are obtained by introducing a 
Lorentzian broadening of 0.4i into (5-functions (vertical bars) 
in the ED results. 
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III. RESULTS AND DISCUSSIONS 

As mentioned in Sec. I, the nature of the ground state 
in the J-J' Heisenberg model changes from the Neel 
to collinear states at J'/ J ~0.6.t§&3 By using the ED 
method, we found that the single-hole spectral function, 
also shows different behavior between the two states.E^I 
However, since we are interested in hole dynamics in the 
AF spin system with realistic values of J'/ J, we focus on 
the region J'/J<0.6. 



A. Zero Temperature 

Figure [j] shows A(k, to) of the t-t'-t"-J model with var- 
ious values of </' at T=0. Parameter values are set to 
be J=0.4t, t'=-0.34t and t"=0.23t as used in Ref. 2. 
Solid lines and vertical bars are the ED results on a 
4x4 cluster. In Fig. 1(a) (J'/J=0), the height of the 
QP peak (w/t~— 1.2) at (ir,0) is small as compared with 

those at (7r/2,7r/2) and (tt/2,0) (u/t 2.2 and -1.7, 

respectively) .□ Since t' and t" terms cause the shift of 
the QP peak at (ir,0) to higher binding energy, the po- 
sition of the QP peak enters into the region where the 
incoherent states exist. As a result, the QP weight is 
strongly suppressed. 

With increasing J' the ED results show the decrease 
of the QP weight for all the momenta. In particular, the 
remarkable decrease is found at (0,0) and (7r/2,0). The 
spectrum near the QP energy at (7r/2,0) becomes broad 
as clearly seen in Figs. 1(c) and (d). Such a broad spi 
trum has been observed at (7r/2,0) in ARPES data! 
Since the QP peak without J' is as sharp as that at 
(7r/2,7r/2) (see Fig. 1(a)), the observed broad spectrum 
is considered to be a manifestation of the effect of frus- 
tration due to </'. 

In order to clarify the effect of J' on the spectral func- 
tion, we also use the SCBA method for the calculation of 
A(k,oj). In Fig. 1, the SCBA results are compared with 
the ED ones to see the validity of the approximation. The 
SCBA data (dashed lines) are plotted for momenta inside 
the magnetic first Brillouin zone. The spectra obtained 
by the SCBA method are in good agreement with those 
by the ED one, except for J' /J— 0.5 where the Neel state 
is expected to be less stable in the ground state. Encour- 
aged by the good agreement between the ED and SCBA 
results at J'/J<0.3, we extend the size of cluster for the 
SCBA method from 4x4 to 16x16. Figure 2 shows the 
SCBA results at various momenta for the 16x16 cluster. 
Shown in the figure are A(k, w), the imaginary part of 
the self-energy ImE (k, uj) and the QP weight Z (k, £qp) 
defined as 



(a) 



Z(k,e QP )= 1 



<9Re£ (k,w) 



where £qp is the QP energy at k. 
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FIG. 2. Single-hole spectral function A(k, tj) (solid line), 
the imaginary part of the self-energy ImE (k, uS) (dashed line) 
and the QP weight Z (k, £qp) (vertical line) of the t-t'-t"-J 
model with J' at T=0 obtained by the SCBA method on a 
16x16 cluster. t'=-0.34t, t"=0.23t and J=0.4f. J'/ J is 
(a) and 0.3 (b). r]/t=0A for A(k,u) and ImE (k, id), and 
?7/i=0.005 for Z(k,e Q p). In (b), Z at (0,0) and (tt/4,0) is 
zero. 

With increasing J'/J from (Fig. 2(a)) to 0.3 
(Fig. 2(b)), the spectral intensity around the QP energy 
becomes broader when the momentum k is close to (0,0). 
Since the life time of the QP expressed by ImE (k, ui) 
shows little change, the broadness comes from the de- 
crease of the QP weight Z. For k=(0,0) and (tt/4,0), Z 
is zero at J'/J=0.3. On the other hand, for large mo- 
menta such as (37r/4,7r/4) and (7r,0), the change of Z is 
very small. In the intermediate values of k, the ratio of 
Z for J' /J— 0.3 to that for J'/J=0 decreases as k ap- 
proaches (0,0) point. The maximum value of the ratio 
appears at (7r/2,7r/2). Similar k dependent behaviors of 
Z are obtained when the value of J in the t-J model is 
reduced. Ej The similarity between the effects due to J in 
the t-J model and J' in our model is also seen in the de- 
pendence of the QP band width W, i.e. with increasing 
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J', W decreases monotonically, which resembles a well- 
known behavior of Wot J in the t-J modelEj Therefore, 
we can say that the changes of A(k, u>) due to the increase 
of J' are roughly similar to those due to the decrease of 
J in the t-J model. This is consistent with an intuitive 
picture that J' frustrate the Neel state as if the net value 
of J is reduced. 
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FIG. 3. Single-hole spectral function A{U.,uj) of the 
t-t'-t" -J model with J' at various temperatures. t'=— 0.34t, 
t"=0.23f and J=0.4t. J'/J is (a) and 0.3 (b). The val- 
ues of temperature are T=0 (solid line), J/2 (dashed line), J 
(dotted line) and t (dot-dashed line). Calculations are per- 
formed on the 4x4 cluster by the ED method. A Lorentzian 
broadening of 0.4t is introduced into the 5-functions. 



B. Finite Temperatures 

We next present the temperature dependence of 
A(k, lo) on a 4x4 lattice. Since the energy separation 
between the ground and first excited states in the J- 



J' Heisenberg model on the lattice is roughly J/2, reli- 
able results are obtained only in the region of T > J/2. 
Although the temperature one can treat is rather high, 
the effect of the change of spin correlation length £ on 
the spectrum may be seen in the calculation because the 
correlation length £ for the temperatures above 
is less than the maximum length of the 4x4 lattice 
(~ 2 v / 2a).E3 

Figure 3 shows A(k, u) of the t-t'-t" -J model with 
J'/J=0 and 0.3 at various temperatures obtained by 
using the ED method. The QP spectral intensity at 
(7t/2,7t/2) decreases with increasing temperature. This 
is not due to the decrease ofpfche QP weight but due to 
the decrease of the life time.Ba We find that, when the 
temperature changes from T=0 to J/2, the intensity in 
Fig. 3(b) (J'/J=0.3) is largely suppressed as compared 
with that in Fig. 3(a). This simply comes from the fact 
that the frustration J' disturbs the AF spin background 
which is responsible for the formation of QP and makes 
the background sensitive to a slight change in temper- 
ature. Such a sensitive behavior of the QP peak may 
be one of the origins of the observed strong temperature 
dependence of ARPES dataM 

In Fig. 3(b), we find that the dependence of the spec- 
tra near the QP energy on temperature at (ir/2,ir/2) is 
remarkably different from that at (7r/2,0). At (ir/2,w/2) 
the spectral weight around u>^—t slightly increases with 
increasing temperature in spite of the reduction of the 
QP intensity. On the contrary, the weight in the same 
energy region at (7r/2,0) is uniformly suppressed with 
increasing temperature. Such a different behavior has 
been observed in ARPES experiment on S^CuC^C^.Iij. 
Since the behavior can not be clearly seen without J' (see 
Fig. 3(a)), the agreement between theory and experiment 
indicates the importance of J' on the electronic states of 
insulating cuprates. 

With increasing temperature from T=0 to J, the QP 
spectral intensity in Figs. 3(a) and (b) rapidly decreases 
accompanied by a slight shift to lower-energy side.LZrE£l 
At T=J, the QP peaks almost vanish because .the corre- 
lation length £ is less than one-lattice spacingO Similar 
suppression of the weight as the temperature increases 
from T=0 to J is seen in the other peaks at u) ~0, t, 2t 
and it for (0,0), {n,ir), (tt,tt/2) and (tt,0), respectively. 
The weight at these peaks is thus sensitive to the nature 
of the spin background. At T—t, the spectrum has a 
broad maximum at a certain energy in each momentum. 
This behavior is independent of the value of J'. In such a 
high temperature region, the hole can pmpagate without 
being disturbed by the spin background. El This produces 
a dispersive structure whose energy is determined by a 
tight-binding band due to the free hole propagation. The 
energy band is given by 

-Efrcc (k) = 2t (cos k x + cos k y ) + 4t' cos k x cos k y (17) 
+2t" (cos 2k x + cos 2k y ) . 

In Fig. 3, the energy obtained from the band is indicated 
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by arrow at each k point. With further increasing tem- 
perature, the maximum of the spectrum is expected to 
shift to the position indicated by the arrow. 



IV. CONCLUSIONS 

We have examined the effect of the magnetic frustra- 
tion J' on the single-hole spectral function in the t-t'-t"- 
J model, by using the ED and SCBA methods. At zero 
temperature, the frustration suppresses the QP weight 
and makes the spectrum broad for small momentum, 
whereas the peak at (7r/2,7r/2) remains sharp. Our re- 
sults explain the experimental data that the (7r/2,0) spec- 
trum in Sr2Cu02Cl2 is broader than the (7r/2,7r/2) one. 
We have also evaluated the temperature dependence of 
the single-hole spectral function in the t-t'-t"-J model 
with J' by using the ED method. The lineshape between 
(7r/2,0) and (7r/2,7r/2) spectra behaves differently with 
increasing temperature. This difference is also consistent 
with recent experimental data. These findings indicate 
the importance of the magnetic frustration on the elec- 
tronic states of the insulating cuprates. 
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